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(54) Process for determining characteristics of suspended particles 



(57) The invention provides an improved process for 
determining the characteristics of dispersed particles. 
The process involves directing waves (acoustic or light) 
into a dispersion, and measuring the attenuation of the 
waves for particular frequencies to provide an attenua- 
tion spectrum. The measured attenuation spectrum is 



then compared to a set of theory-based calculated at- 
tenuation spectra to determine the particle size distribu- 
tion corresponding to the measured attenuation spec- 
trum. Unlike previous processes, the particle size distri- 
bution is capable of being accurately determined by a 
single inversion algorithm, thus eliminating burdensome 
steps of previous methods. 
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Description 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

[0001] The invention relates to characterization of dispersed particles, e.g., determination of particle size distribution. 
Discussion of the Related Art 

10 

[0002] Processes such as chemical-mechanical planarization (CMP) and sol-gel fabrication of silica bodies utilize 
dispersions of particles in a liquid medium. (For a discussion of CMP and sol-gel processes, see, respectively, P. Van 
Zant, Microchip Fabrication, 3rd Ed., McGraw-Hill, 300 (1997), and U.S. Patent No. 5,240,488.) For such processes, 
it is desirable to know the specific properties of the particles. For example, larger particles in CMP often lead to scratches 
f 5 and variations in the median particle size and the size distribution sometimes lead to inconsistent polishing. Similarly, 
the properties of sol-gel derived bodies depend in large part on the characteristics of the particles in the dispersion. 
For these reasons, among others, various methods have been developed to determine the characteristics of particles 
in such dispersions. 

[0003] One such characteristic is the particle size distribution (PSD). See, for example, U.S. Patents Nos. 4,706,509, 
20 5,121,629, and 5,569,844 (the '509, '629, and '844 patents, respectively). The '509 patent relates to a method for 
ultrasonically measuring solids concentration and particle size distribution in a dispersion. Ultrasonic waves at a variety 
of frequencies are directed into the dispersion, and the attenuation at the frequencies is measured. A dimensional 
spectrum (across the range of particle dimensions) is divided into dimensional intervals, and a system of linear equa- 
tions is developed to represent the concentration of particles in each dimension interval. The system of equations is 
25 then solved to determine the PSD. However, the process of the '509 patent is limited to waves having a wavelength 
smaller than the largest particles. The process is therefore not useful for relatively fine particles, e.g., tens to hundreds 
of nanometers, since acoustic waves cannot presently be driven at the necessary frequencies. 
[0004] The process of the '844 patent involves measuring the attenuation of both ultrasonic waves and electromag- 
netic radiation to determine particle size distribution. Specifically, ultrasonic velocity and ultrasonic attenuation are 
30 combined with the density, as determined from the electromagnetic radiation attenuation, to calculate the PSD. In 
addition to the problems inherent in generating x-rays or gamma rays, however, particle sizes of about 10 to 15 urn 
appear to be the lower limit for the process. 

[0005] In the '629 patent, the disclosure of which is hereby incorporated by reference, ultrasonic waves at a variety 
of selected wavelengths are passed through a dispersion, and the attenuation at each frequency is measured to derive 
35 a measured attenuation spectrum over those frequencies. Separately, based on a theoretical model, a set of attenuation 
spectra are calculated for a variety of PSDs, and the calculated spectra are then compared to the measured spectrum 
to formulate a preliminary approximation of the PSD of the dispersion. Further calculations must be performed, starting 
from this approximation, to more accurately determine the PSD. 

[0006] While current techniques for characterizing dispersed particles, e.g., determining PSD, are adequate for many 
40 purposes, improvements that provide more accurate, consistent, and reliable results are desired. 

Summary of the Invention 

[0007] The invention provides an improved process for determining the characteristics of dispersed particles. (As 
45 used herein, the term particles includes solids, liquids, or gases dispersed in a continuous medium.) The invention 
involves directing waves (acoustic or light) into a dispersion, and measuring the attenuation of the waves for particular 
frequencies to provide an attenuation spectrum. (As used herein, waves includes, but is not limited to, continuous 
waves and pulsed measurments.) The measured attenuation spectrum is then compared to a set of theory-based 
calculated attenuation spectra to determine the particle size distribution corresponding to the measured attenuation 
50 spectrum. Unlike previous processes, the particle size distribution is capable of being accurately determined by a single 
inversion algorithm, thus eliminating burdensome steps such as those of the '629 patent, e.g. deriving a preliminary 
approximation of PSD, and going forward from the approximation through several inversion techniques. (Inversion 
techniques, in the present context, involve taking a set of known particle size distributions, determining the attenuation 
spectrum that each PSD would theoretically produce, and comparing a set of such theory-based spectra to the actual, 
55 measured spectrum to find the actual PSD.) By avoiding steps such as an initial approximation and a multi-step inver- 
sion technique, the invention produces more accurate, consistent, and reliable results. 

[0008] In particular, while the process of the '629 patent provides useful results for some applications and some 
particle size distributions, the process has several potential drawbacks. The '629 patent starts with a well-known model 
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for acoustic attenuation, which was set forth in Epstein and Carhart, J. Acoust. Soc. Amer., Vol. 25, 1953, 553, and 
revised by Allegra and Hawley, J. Acoust. Soc. Amer., Vol. 51, 1972, 1545, the disclosures of which are hereby incor- 
porated by reference. From this model (referred to herein as the Allegra-Hawley model), a set of attenuation spectra 
are calculated for a set of known particle size distributions. When the actual attenuation spectra of a dispersion is 

5 measured, the '629 patent takes an initial guess of the PSD based on an initial comparison of the calculated spectra 
and measured spectrum (see Col. 1 5, lines 30-39). The quadratic discriminator (the first of three inversion techniques) 
then attempts to refine this initial guess by computing a limited number of attenuation spectra similar to the initial guess. 
Following the quadratic discriminator, the'629 patent uses the Powell discriminator, which selects a specific functional 
form with a limited number of parameters, in an effort to speed the process of finding a match to the measured atten- 

10 uation spectrum. However, this inversion technique may generate an inaccurate response when the selected functional 
form is not flexible enough to describe the shape of the true PSD. Finally, where the initial guess in the '629 patent 
indicates a narrow PSD, the impulse discriminator inversion technique is used, this third technique using linear equa- 
tions to determine a PSD with effectively no spread in particle size. Unfortunately, the applicability of this third technique 
is limited. Specifically, because the impulse discriminator appears to be tailored toward specific systems, it has the 

15 potential to generate inaccurate responses when a variety of systems are tested. Moreover, the process of the '629 
patent is tailored to acoustics, and would not be useful for light attenuation. 

[0009] The invention overcomes these drawbacks of the '629 patent by application of a faster, more straightforward 
inversion technique. Specifically, the invention uses a single inversion algorithm that requires no prior knowledge (or 
guess) of the particle size distribution. Moreover, the single inversion algorithm is independent of the physical meas- 
20 urements being made relating to particle characteristics, i.e., the inversion technique is applicable to any model for 
predicting attenuation based on PSD, such as an acoustic attenuation model or a turbidity, i.e,, light attenuation, model. 
In one embodiment, an improved Allegra-Hawley model is used to calculate the set of theory-based acoustic attenuation 
spectra. 

[0010] The single inversion algorithm involves parameterizing a particle size distribution, i.e., mathematically de- 
25 scribing an abstract form, such as a PSD based on a selected number of parameters. The parameters are then varied 
to obtain a "best" agreement between a measured spectrum and theoretical spectra, this best agreement providing 
the PSD. The best agreement is defined by minimization of an error function that is selected to not only incorporate 
the differences between measurement and theory, but to also weight each measured data point according to the con- 
fidence in that data. Unlike the *629 patent, no initial guess of the PSD is needed. 
30 [0011] The invention therefore provides an improved process for determining the characteristics, particularly particle 
size distribution, of a particle dispersion, including a dispersion containing nanometer-size particles. 

Brief Description of the Drawings 

35 [0012] Fig. 1 schematically illustrates an apparatus useful in the invention. 

[0013] Fig. 2 illustrates the improvements made in the Allegra-Hawley model, according to the invention. 

[0014] Figs. 3A and 3B show particle size distributions generated according to the invention from simulated acoustic 

attenuation spectra. 

[0015] Figs. 4A and 4B show, respectively, measured acoustic attenuation data and a particle size distribution gen- 
40 erated from the data by the process of the invention. 

[0016] Figs. 5A-5C show, respectively, measured acoustic attenuation data, a particle size distribution generated 
from the data by a previous method, and a particle size distribution generated from the data by the process of the 
invention. 

[0017] Figs. 6A and 6B show, respectively, raw turbidity data and a particle size distribution generated from the data 
45 by the process of the invention. 

Detailed Description of the Drawings 

[0018] Acoustic attenuation techniques for characterizing particles in dispersion involve the interaction of applied 
50 sound waves with the dispersed particles. As a sound wave travels through a dispersion, the wave loses acoustic 
■ energy by various scattering mechanisms. The total energy loss is given by the attenuation coefficient, a (dB/inch): 

a=Jln^ (1) 
55 L '1 

where L is the path length across the dispersion, l 0 is the initial intensity of the sound wave, and. I., is the intensity of 
the sound wave after passing through the dispersion. Measurement of the attenuation coefficient at different frequencies 
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of the sound wave leads to an acoustic attenuation spectrum. The above-referenced models by Epstein and Carhart, 
and by Allegra and Hawley make it possible to predict the attenuation spectrum for particles of a given size distribution 
and concentration. The models require knowledge of several physical properties of the particles and the liquid medium, 
including density, thermal expansion coefficient, thermal conductivity, heat capacity, viscosity, and shear rigidity. It is 

5 then possible to construct a 3-D matrix that relates attenuation, frequency, and particle size. 

[0019] Acoustic attenuation in a particle dispersion is measured by use of an apparatus such as the Ultrasizer™, 
made by Malvern Instruments, Ltd., Worcestershire, United Kingdom (who is also the assignee of U.S. Patent No. 
5,121,629, discussed above). A schematic of the chamber 10 of such an apparatus is shown in Fig. 1. Two pairs of 
broadband transducers 12, 14, 16, 18, are in contact with a sample located in a tank 20, typically formed of stainless 

10 steel. Advantageously, the transducers are designed to transmit a wave as close to a plane wave as possible. The 
chamber 10 typically contains a temperature sensor, heater, and a propeller (none of which are shown). It is possible 
to remove the sample by use of a drain valve at the base of the tank 20. Typically, one pair of transducers 12,14 covers 
a lower frequency range, e.g., 1 to 20 MHz, and the other pair 16, 18 a higher range, e.g., 15 to 200 MHz. Generally, 
the transmitting transducers 12, 16 are capable of being moved to different positions in the chamber, but the receiving 

15 transducers 14, 18 are fixed. The spacing between the transmitting transducers 12, 16 and the receiving transducers 
14, 18 is controlled by a stepper motor. The transmitting transducers 12, 16 direct a series of planar sound waves, 
generally at preselected frequencies, across the sample and the receiving transducers detect the attenuated waves. 
Computer and signal processing equipment, as apparent to one skilled in the art, is used to control the process of 
attaining the attenuation spectrum, monitor the detected sound waves, and manipulate the signals generated therefrom. 

20 As commercially available, the Ultrasizer™ provides both raw attenuation data as well as predicted particle size dis- 
tributions. Other apparatus that measure attenuation, such as discussed in references cited in the Background section, 
are equally useful in the invention. Other commercially-available apparatus capable of being used in the invention 
include the CMP Slurry Monitor, made by Colloidal Dynamics, Inc. of Warwick, Rhode Island, and the AcoustoPhor 
device, made by PenKem. It is also possible to use devices that utilize Fourier transform techniques to measure at- 

25 tenuation spectra. Advantageously, the apparatus is capable of being placed in-line in a production process, such that 
real-time measurements are provided. 

[0020] Optical attenuation, or turbidity, measures the loss of intensity of light as it travels a fixed distance through a 
suspension. A turbidity spectrum is obtained by measuring attenuation at several wavelengths (frequencies) of light. 
The Mie theory describes the basic interaction between light and a dispersion of spherical particles. (See G. Mie. Ann. 
30 Physik , Vol. 25, 377 (1908). See also M. Kerker, The Scattering of Light and other Electromagnetic Radiation, Academic 
Press (1969).) The loss of intensity is directly related to the difference in refractive indices of the particle and the liquid 
medium, as well as to the particle size distribution. The attenuation mechanisms for light and sound differ greatly, 
making independent measurements by these techniques complementary. 

[0021] The measurement of optical attenuation is able to be performed using a conventional ultraviolet-visible region 
35 (UV-Vis) spectrophotometer. Commercial units are available from several manufacturers such as Perkin-Elmer Corp., 
Spectral Instruments Inc., and Ocean Optics Inc. A basic device includes a broadband light source, a sample cell, a 
UV-Vis monochromator and a photodetector. Light from the broadband source is collimated and directed through the 
sample cell. The light passing directly through the sample cell arrives at the monochromator, which directs a finite width 
of frequencies to the photodetector. By adjusting the monochromator, the light transmitted by the sample can be meas- 
40 ured as a function of frequency, the log of the ratio of this spectrum to a standard reference is the UV-Vis attenuation 
or absorption spectra. 

[0022] A variety of dispersions (solid, liquid, or gas particles in a continuous medium) are capable of being measured 
according to the invention. For example, it is possible to measure solid-in-liquid dispersions, such as silica, alumina, 
titania, or other oxides, in aqueous or non-aqueous media. It is also possible to measure emulsions (i.e., liquid particles 
45 in a continuous liquid medium), as well as liquid-in-gas and gas-in-liquid dispersions. The particle concentrations and 
particle diameters for which the invention will be functional are capable of varying widely, as long as a useful attenuation 
spectrum is obtained. Typical concentrations of particles that provide useful results range from 0.1 to 60 vol.%, more 
typically about 1 to 60 vol.%, and typical particle diameters range from about 1 nm to about 1000 |im, more typically 
about 10 nm to about 10 nm. 

so [0023] The invention involves directing waves (acoustic or light) into a dispersion, and measuring the attenuation of 
the waves for each particular frequency to provide an attenuation spectrum. The measured attenuation spectrum is 
then compared to a set of theory-based calculated attenuation spectra to determine the particle size distribution cor- 
responding to the measured attenuation spectrum. The comparison is performed by a single inversion algorithm which 
requires no prior knowledge (or guess) of the particle size distribution, and which is independent of the nature of the 

55 physical measurements being made, e.g., acoustic or turbidity. The algorithm involves parameterizing a particle size 
distribution, and then varying the parameters to obtain a "best" agreement between measurement and theory, as de- 
fined by minimization of an error function that is selected to' incorporate the differences between measurement and 
theory, and also weight each measured data point according to the confidence in that data. The minimization of this 
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error function with respect to variations in PSD is capable of being performed through a conjugate gradient minimization 
procedure. The details of the process are as follows. 

Inversion Technique 

5 

(1) Parameterizing a particle size distribution 

[0024] As discussed above, parameterizing involves mathematically describing an abstract form, such as a PSD, 
based on a selected number of parameters. In the invention, the particle size distribution is initially viewed as a corn- 
to pletely general form, <j> (x), which represents particles of varying size x. No assumptions or guesses as to the form of 
the PSD are made. 

(2) Defining the Error 

15 [0025] The error function, E, numerically describes the agreement between the theoretical and experimental spectra. 

E = \dfW(f)[P(f)-0(f)] 2 (2) 

20 Here, f represents frequency, W(f) is a weighting function, P(f) is the predicted spectrum, and O(0 is the observed 
experimental spectrum. This form is similar to a standard least-squares formulation where the square ensures non- 
negativity. Defining E as an integral over continuous functions rather than referring directly to the discrete frequencies 
of the experimental spectrum ensures that E's value remains independent of the details of the experiment (i.e. the 
number and location of the discrete frequency measurements). 

25 [0026] The inclusion of the weighting function, W(f), accounts for frequency dependent experimental uncertainty. For 
example, higher frequency data tends to include more experimental uncertainty (e.g., noise) than low frequency data. 
In such a case, giving the low frequency data more weight would emphasize the importance of differences in the low 
frequency data. Given a suitable standard, it is also possible to include instrumental variations within this weighting 
factor. Although the optimal W(f) depends upon experimental conditions and the exact nature of the measurements, 

30 a nearly optimal form is obtained by weighting the data based upon the variation of the experimental data over several 
repeated measurements, 

W{f } )=\ (3) 

35 °j 

where cj represents such statistical variations, with j indicating particles in the dispersion. 

[0027] The predicted spectrum, P(f), is given by a linear combination of the attenuation from each particle within the 
distribution. 



40 



P(f) = \dRA(f,RMR) (4) 



Although this predicted spectrum is an approximation (owing to nonlinearities induced by multiple scattering phenom- 
45 ena), it provides desirable results in the invention. (The fact that this spectrum is an approximation, however, does not 
constitute an initial guess of PSD.) Within this approximation, the function A(f,R) describes, in its entirety, the physical 
attenuation model. The function represents the attenuation at frequency f for a single particle of size R. By changing 
the form of this function, we can use this inversion method for both optical and acoustic attenuation. 
[0028] Further, the treatment of models which are non-linear in §(R), e.g., those for treating multiple scattering phe- 
50 nomena, may be accommodated by adding terms to Equation 4 representing higher powers of §(R), without further 
changes to the invention. 

[0029] Initially, as mentioned above, the particle size distribution is viewed as a completely general form, $(/?), unlike 
previous methods. A physically reasonable function form for §(R), i.e., an initial approximation of the PSD, could be 
selected as in such previous methods. However the invention avoids such an approximation in favor of a direct dis- 
55 cretization of the above equations (discretization is the representation of a continuous function through the values of 
that function at discrete points). Therefore, <)>(/?) is represented as its magnitude at discrete values ty, and thus the 
representation of ${R) becomes increasingly complete as the number of discrete particle sizes treated increases. 
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(3) Computing the Error Function 

[0030] Because the measurement of attenuation is performed at discrete frequencies, the computational treatment 
of this problem requires the discretization of the error function. In addition, because the attenuation data tends to be 
5 logarithmically spaced, the discretization is also desirably logarithmic. In this case the integrals (Eqns. 2 & 4) become 
summations, 

£ = Z(Alg/)/^[^-0,] 2 (5) 



/>=£(Alg/?)K,V> (6) 

where / represents frequency and j represents particle sizes. 

Owing to the exponential ranges of frequencies and particle sizes, the terms of these summations may vary over 
20 several orders of magnitude. This means the accurate computation of E (and P,) requires careful attention to round- 
off accumulation in the summations. 

[0031] The representation of the theoretical model is contained within the matrix^, which is optionally pre-computed 
for particular physical models, e.g. the acoustic attenuation of silica in water. This pre-computation allows the time- 
consuming portion of the process (i.e. model evaluation) to be performed independent of the measurement and inver- 
ts sion process. Since the experimental measurement defines the discrete frequencies of interest, the/^y used for a given 
inversion is often obtained through interpolation of a dense, pre-computed A t y 

[0032] For a given experimental measurement, the discrete frequency values remain fixed throughout the inversion 
process. Thus, for computational speed, it is possible for the frequency summations to be pre-computed, leading to a 
more quickly computed error. 

30 

E^2( A ^ R ) 2]B ^A + y £{MzR)C i <t> j +E 0 (7) 

jk j 

35 



B ik = X(Aig/)/^(Ai g /?) 2 AA W 

i 

40 



Cj = ^(Alg/^Alg/?)^ (9) 

45 i 



Although these equations are formally identical to Equation 5, numerical round-off leads to computational differences. 
The invention uses both forms of the error expression to balance speed and accuracy. 

55 

(4) Minimizing the Error 

[0033] The determination of the particle size distribution requires minimization of the error function with respect to 
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variations in the ty. Physically, the particle size distribution must be non-negative for all particle sizes. This introduces 
constraints on the possible variations in the fy. To remove this constraint, the invention introduces another set of var- 
iables Cj such that 

$j = (Cj) 2 . (11) 

This ensures the are non-negative. Thus, instead of requiring a difficult constrained minimization onty, it is possible 
to perform a simple unconstrained minimization on the independent, unconstrained Cy. 

10 [0034] The minimization begins with an unbiased starting distribution where all Cy are zero. To speed convergence 
from this starting point, several sub-minimizations are performed. The sub-minimizations correspond to coarse discre- 
tizations of the §(R). For example, it is possible for the first sub-minimization to divide the range of possible particle 
sizes into 50 discrete values. Upon reaching a suitable minimal error for this sub-minimization (defined when the gra- 
dient of the error is below a chosen tolerance), the number of discrete values is increased (e.g. to 100 particle sizes). 

15 The information from a previous sub-minimization is used as an initial starting point for subsequent minimizations. This 
process continues until the final discretization level is achieved. This allows the code to quickly select a region of 
particle sizes in the early sub-minimizations, which is further refined in future minimizations. 

For this minimization, the invention employs a conjugate gradient minimization using exact line searches. The search 
directions are accurately computed using the gradient of Equation 5. The line search is quickly processed using the 
20 gradient of Equation 7. It is possible that scaling of the gradient, or other pre-conditioning, will increase the speed of 
the minimization. 

Attenuation Model 

25 [0035] As mentioned above, the inversion technique of the invention is independent of the physical measurements 
being made to characterize the particles. Thus, there are at least two models capable of being combined with the 
inversion technique to determine particle size distributions. The first model is that initially developed by Epstein and 
Carhart, referenced above, relating to acoustic attenuation of water droplets in fog, and then further developed by 
Allegra and Hawley in the context of solid spherical particles in continuous liquid media. However, to improve the results 

30 obtained by the invention, several modifications have been made to this standard Allegra-Hawley formulation (referred 
to as improved Allegra-Hawley). These modifications concern what is believed to be an unjustified approximation for 
the treatment of solid particles, as well as several smaller changes to the form of the equations necessary to obtain a 
numerically stable and thermodynamically consistent computational method. Although each of these changes is easily 
understood by one skilled in the art of numerical computing, it is important to note that, without such careful attention 

35 to numerical stability and round-off, the invention would not perform as well. 

[0036] Many of the modifications to Allegra-Hawley are computational rather than mathematical. The large ranges 
of particle size, frequency, and attenuation apparent for the acoustic measurements require careful attention to numer- 
ical stability. The basic computational task of obtaining a theoretical attenuation spectrum requires the formation of a 
set of coefficients and the solution of a set of coupled linear equations. The computation of these coefficients requires 

40 the evaluation of several auxiliary functions of complex argument. Care must be taken that these functions are reliably 
computed over the range of arguments necessary. Also, the standard numerical tasks of minimizing round-off and 
ensuring the linear equation solver is robust and stable are important for correct computation of theoretical spectra. 
This attention to detail allows the invention to compute numerically correct attenuation spectra for, at a minimum, 
particles ranging in radius from 5 nm to 10 ^im and for frequencies from 100 kHz to 200 MHz. The approximations 

45 inherent in the physical model, rather than numerical issues, tend to increase the potential for errors as these ranges 
are extended. 

[0037] Thermodynamics provides exact relationships between the physical properties of a material system. Ther- 
modynamic consistency refers to methods which maintain these exact relationships. For computational methods, ther- 
modynamic consistency is manifested as stability with respect to minor changes in input physical properties. With the 

50 Allegra-Hawley formulation, researchers have found the computed attenuation spectrum to be very sensitive to small 
changes in the input quantities. By re-deriving the equations using quantities that remain thermodynamically consistent, 
our computational method does not show these sensitivities to such input properties. One important example of ob- 
taining thermodynamic consistency concerns the relationship between the specific heat ratio, the isothermal compress- 
ibility, and the speed of sound. The specific heat ratio, y. is determined by the speed of sound q, the coefficient of 

55 thermal expansion p, the temperature T, and the specific heat at constant pressure C p . 
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Y =1 + 2JJ (12) 
5 The isothermal compressibility, k, is determined by the speed of sound, the specific heat ratio, and the density p. 

* = (13) 
Q P 

10 

These two equations provide an exact relationship between y , q, and k. It is possible, however, that independent 
experimental determinations of y, q, and k will not satisfy these equations exactly Small deviations away from these 
exact relationships become important since these quantities enter the final equations as the difference of two large 
numbers. Thermodynamic consistency is achieved by using numerically stable intermediates satisfying the constraints 

15 inherent in the exact relationships. 

[0038] As an illustration, Table 1 shows the average changes in the theoretical acoustic attenuation (for R = 5 nm to 
10 ^im; and for f = 1 MHz to 150 MHz) given a ± 10% change in each physical property. These resulting changes in 
attenuation are small relative to the range of measured attenuations, andthis table therefore illustrates the stability of 
the invention with respect to changes in input properties, where this stability is due to the re-derivation of the model's 

20 equations. As expected, the table also shows that the viscosity, density, and speed of sound are the most important 
physical characteristics of the system for predicting attenuation. The stability in attenuation provided by the invention's 
use of an improved Allegra-Hawley leads to a more accurate determination of particle size distribution. 



Table 1 



Property i 


Atten. (+10% change in property) 


Change (-10% change in property) 


Particle: 


Specific Heat 


0.1 % 


0.1 % 




Density 


15.1 % 


18.0% 




Speed of Sound 


0.3 % 


0.4 % 




Thermal Dilation 


-0.0 % 


-0.0 % 




Lame Coefficient 


0.1 % 


0.0 % 




Attenuation (pure) 


-0.0 % 


-0.0 % 


Liquid: 


Specific Heat 


0.1 % 


0.1 % 




Density 


19.2% 


18.4% 




Speed of Sound 


12.0% 


10.3% 




Thermal Dilation 


0.1 % 


0.1 % 




Viscosity 


4.3 % 


5.0 % 




Attenuation (pure) 


0.8 % 


0.7 % 



[0039] The improved Allegra-Hawley is also based on the following: The standard Allegra-Hawley formulation defines 
the compressional wavevector k c as 

k e = J + ia L (14) 

where co is the applied frequency, q is the speed of sound and a L is the acoustic attenuation for the pure substance. 
This expression is obtained by neglecting second order terms when decoupling the compressional and thermal wave 
equations. While this approximation is justified in the case of liquids or gases (as in Epstein & Carhart's work), for solids 
this approximation incorrectly computes the compressional wavevector. Fortunately, this approximation is easily re- 
moved from the standard formulation by computing the compressional and thermal wavevectors exactly, including 
second and higher terms. Figure 2 displays the effect (in percentage change in attenuation) of removing this approx- 
imation from the computed attenuation spectra for three particle sizes (radii of 25 nm, 165 nm, and 500 nm). The effect 
is larger for larger particles at higher frequencies. Because the change in attenuation varies with frequency, this error 
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would lead to improper inversions. 

[0040] The second model useful in the invention relates to turbidity, which is the attenuation of light, and was initially 
presented by G. Mie, Ann. Physik , Vol. 25, 377 (1908). See also M. Kerker, The Scattering of Light and other Electro- 
magnetic Radiation, Academic Press (1 969). This model is capable of being successfully used in the invention without 
5 any modifications. 

[0041] The invention will be further clarified by the following examples, which are intended to be exemplary. 
EXAMPLE 1 

10 [0042] The inversion technique and the model were first tested in the absence of experimental error using a set of 
simulated attenuation spectra. The simulated spectra were obtained from the acoustic attenuation model for two log- 
normal particle size distributions fti=25 nm (radius), o=0.005; n=500 nm (radius), o=0.1) The inversion process, when 
applied to the simulated spectra, generates a particle size distribution which may be directly compared with the assumed 
distributions. Fig. 3A shows the assumed and computed distributions for both the 25 nm radius and the 500 nm radius 

15 tests. The invention was able to reproduce both the median radius and width of the assumed distributions accurately 
over a broad particle size range. 

[0043] Simulated spectra were also used to test the sensitivity of the invention toward the presence of more than 
one population of particle size. Fig. 3B shows a set of particle size distributions obtained by inverting simulated spectra 
for three different mixtures of 25 nm radius and 165 nm radius particles (again superimposed on the assumed distri- 
20 butions). The ratios in the Figure refer to the volume percentages of 25 nm radius and 1 65 nm radius particles, respec- 
tively, in each mixture. The two populations of particle sizes were distinguished clearly by the invention, and the relative 
proportions and widths of the distributions were reproduced with reasonable accuracy. 

EXAMPLE 2 

[0044] An acoustic attenuation spectrum was measured for a 1 vol% suspension of silica particles (MP4030, Nissan 
Chemical Corp., Tokyo) using an Ultrasizer™ (Malvern Instruments Ltd., Worcestershire, U.K.). The tested Nissan 
particles are often used as a reference material for sizing instruments, having a nominal diameter of 300 nm and a 
relatively narrow size distribution. The sample was stirred at 150 rpm with a stainless steel propeller, and the temper- 
30 ature was 25°C. The automatic optimal strategy feature of the Ultrasizer™ was used to determine the most effective 
frequency range and transducer spacings for the attenuation measurements. Fig. 4A shows the measured attenuation 
spectrum. The resulting particle size distribution determined by the inversion technique of the invention is shown in 
Fig. 4B. The particle radius (-0.165 |im) was accurate, and the spread of the distribution was narrow, as expected. 

35 EXAMPLE 3 

[0045] Acoustic attenuation measurements were made on commercially available chemical mechanical polishing 
(CMP) slurries. Klebosol 30N50 and Klebosol 30N25 (Solution Technology, Inc., Monroe, NC) are 30 wt% silica slurries 
in ammonia solution, with nominal particle diameters of 50 nm and 25 nm respectively. The 30N50 and 30N25 slurries 
40 ' were each diluted to 5 vol% in water and measured separately in an Ultrasizer™, under the experimental conditions 
described in Example 2. 

[0046] The acoustic attenuation spectra from the Ultrasizer™ for the two samples are plotted in Fig. 5A. The repro- 
ducibility of these spectra was high, with a variance of less than 5% over 60 measurements on each sample. The 
resulting particle size distributions determined by the Ultrasizer™ software are shown in Fig. 5B. These size results 

45 were obtained using the default values for the geometric standard deviation (narrowest allowed deviation = 1 .2, broad- 
est allowed deviation = 3.0) and the default number of iterations (1 5 iterations). The default model for the size distribution 
(a linear combination of two log-normal distributions) was used. Despite a low reported fit error, the 25 nm slurry gave 
misleading size results (of about 500 nm and about 5 ^im) and the estimated concentration (0.8 vol%) was inaccurate. 
The 50 nm sample showed a first mode with mean diameter 41 nm containing 95% of the particle population, with a 

50 second mode of mean diameter 540 nm. 

[0047] Fig. 5C shows the particle size distributions determined by the invention for the measured attenuation spectra 
(Fig. 4A). For both slurries, the main size mode occurred in the vicinity of the nominal particle radius. Narrow peaks 
with small populations were also apparent for both slurries at larger particle sizes. 

55 EXAMPLE 4 

[0048] Turbidity measurements were made using an Ocean Optics S2000 spectrometer (Ocean Optics Inc., Dunedin, 
FL) and a tungsten-halogen light source. The dark current of the photodetector was measured and subtracted from all 



EP 0 989 397 A1 



experimental measurements. The attenuation spectrum of water was used as a reference spectrum for the attenuation 
measurements of the suspensions. Fig. 6A is an example of a measured turbidity spectrum for 50 nm diameter silica 
particles at 1 vol% (Klebosol 30N50, Solution Technology Inc., Monroe, NC). The smooth curve is the theoretical fit to 
the experimental data. Strong absorption (e.g., in the UV region for silica particles) limits the usable wavelength range 
5 in dense suspensions. However, it was found that the 500-900 nm wavelength range was sufficient for the inversion 
process. The particle size distribution determined by the invention for this turbidity spectrum is shown in Fig. 6B. The 
mean particle radius of 21 nm was in accord with the nominal size, and the size spread was narrow. 
[0049] Other embodiments of the invention will be apparent to those skilled in the art from consideration of the spec- 
ification and practice of the invention disclosed herein. 



Claims 



1. A process for determining the characteristics of particles dispersed in a medium, comprising the steps of: 

15 

directing waves into a dispersion; 

measuring the attenuation of the waves for particular frequencies to provide an attenuation spectrum; 
comparing the measured attenuation spectrum to a set of theory-based calculated attenuation spectra to de- 
termine the particle size distribution corresponding to the measured attenuation spectrum, wherein the particle 
20 size distribution is determined in the absence of a step of deriving a preliminary approximation based on a 

comparison between a measured attenuation spectrum and calculated attenuation spectra. 

2. The process of claim 1, wherein the comparison step and the determination of particle size distribution are per- 
formed by a single inversion algorithm. 

25 

3. The process of claim 1 , wherein the comparison step and the determination of particle size distribution are capable 
of being performed with no prior knowledge of particle size distribution. 

4. The process of claim 1 , wherein the waves are selected from acoustic waves and light waves. 

30 

5. The process of claim 4, wherein the waves are acoustic waves and the theory-based attenuation spectra are 
calculated from an improved Allegra-Hawley model. 

6. The process of claim 4, wherein the waves are light waves and the theory-based attenuation spectra are calculated 
35 from the Mie model. 

7. The process of claim 1 , wherein the particles range in diameter from about 1 nm to about 1000 [im. 

8. The process of claim 7, wherein the particles range in diameter from about 10 nm to about 10 |im. 

40 

9. The process of claim 1, wherein the particles are present in an amount ranging from 0.1 to 60 vol.%. 



10. The process of claim 1, wherein the step of comparing comprises: 

45 computing an error function based on comparison between the measured attenuation spectrum and the cal- 

culated attenuation spectra; and 
minimizing the error function. 

11. The process of claim 10, wherein the minimization is performed by conjugate gradient minimization. 

50 

12. The process of claim 10, wherein at least a portion of the computed error function is pre-computed prior to directing 
the waves into the dispersion, based on the physical properties of the dispersion. 



55 
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FIG. 3A 
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FIG. 5A 
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